Statistics of Coulomb blockade peak spacings for a partially open dot 
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We show that randomness of the electron wave functions in a quantum dot contributes to the 
fluctuations of the positions of the conductance peaks. This contribution grows with the conductance 
of the junctions connecting the dot to the leads. It becomes comparable with the fluctuations coming 
from the randomness of the single particle spectrum in the dot while the Coulomb blockade peaks 
are still well-defined. In addition, the fluctuations of the peak spacings are correlated with the 
fluctuations of the conductance peak heights. 
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I. INTRODUCTION 



Quantum dot in the Coulomb blockade regime is con-, 
ventionally described by the constant interaction model.tl 
In this model, the Hamiltonian of the system is repre- 
sented by a sum of two terms. The first one is the elec- 
trostatic charging energy, that does not fluctuate and de- 
pends on the total number of electrons only; the second 
term is the Hamiltonian of free quasiparticles (all inter- 
actions except the charging energy are ignored). For a 
disordered or chaotic quantum dot, free quasiparticles in- 
side the_jdot are described by the random matrix theory 
(RMT).0 The spectrum of quasiparticles is random, with 
the level spacings obeying the Wigner-Dyson statistics. 
The one-particle wave functions are also random, and 
their statistical properties within the RMT theory are 
defined by the Porter-Thomas distribution. 

The predictions of this model may be checked in trans- 
port experiments with a dot weakly coupled by point con- 
tacts to the electron reservoirs (leads). The number of 
electrons on the dot in this setup is nearly quantized, and 
can be controlled by an additional electrode, gate, which 
is capacitively coupled to the dot. The dependence of the 
tunnel conductance through the system, G, on the gate 
voltage Vg exhibits sharp peaks. A Coulomb blockade 
peak corresponds to a point where the ground state of 
the dot is degenerate: the states with n and n + 1 elec- 
trons have the same energy. The (n -I- l)st electron can 
therefore freely tunnel to and from the dot, so G{Vg) has 
a peak at this point. 

For an almost closed dot, the height of the peak is re- 
lated to the values of the wave functions at the points of 
contacts.D The randomness of the wave functions trans- 
lates into the randomness of the peak heights. The exist- 
ing experimental resultsa for the peak conductance dis- 
tribution function are an a satisfactory agreement with 
the theory predictions .0 

In the framework of the constant interaction model, 
the spacings between the conductance peaks contain two 
contributions. The first one is due to the charging energy, 
and does not fluctuate. The second term in the limit of 
weak dot-lead tunneling is proportional to the spacing 
between the discrete energy levels of the RMT Hamilto- 



nian. This term does fluctuate, and obeys the Wigner- 
Dyson statistics. However, the experimentally measured 
distribution of peak spacings is similar to Gaussian ,B~EI in 
apparent disagreement with the predictions of the con- 
stant interaction model. One of the possible explanations 
of this phenomenon is that the electron-electron interac- 
tion, which is not accounted for by RMT, substantially 
affects the dot's spectrum.B Another explanation refers 
to the change in the dot shape with the variation of the 
gate voltage Vg, so the two adjacent peaks in the G{Vg) 
dependence correspoiui to two different realizations of the 
random Hamiltonianjj 

The above-mentioned works concentrated on the al- 
most closed dots, when the tunnel junctions are used 
only as probes, withoutj-affecting the states in the dot. 
In the recent experimentll3 the statistics of peak spacings 
was studied for a partially open quantum dot. In such a 
system, charge quantization is lifted gradually, while the 
conductance of the dot-lead contacts increases. Simulta- 
neously, a crossover occurs between the sharp Coulomb 
blockade peaks and conductance fluctuations of relatively 
small amplitude in an almost open dot. The main part of 
the peak spacing comes from the charging energy. There- 
fore lifting of the charge quantization may yield an ad- 
ditional significant contribution to the peak spacing fluc- 
tuations. This many-body problem was not addressed 
theoretically by now. 

In this paper, a theory of the fluctuations of the peak 
positions in the case of a partially open dot is devel- 
oped. In Sections and III we consider the regime in 
which the conductances of the dot-lead point contacts 
are not negligible, but still considerably smaller than the 
quantum unit Go = e'^/irfi. Under these conditions, the 
conductance peaks remain well-defined. However, their 
positions are already affected by the finite dot-lead cou- 
pling. The randomness of the electron wave functions 
controlling the coupling, contributes to the fiuctuations 
of the peak spacing. We calculate the variance of this 
contribution, and show that it increases with the junc- 
tion conductance. We study the fluctuations of spacings 
between the adjacent, as well as between distant peaks. 
This theory allows us to fojem-a consistent understanding 
of the experimental resultaZlllj for the moderate values of 
junction conductances. 
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In Section [V we consider qualitatively the case of an 



(almost) open dot, when the conductance of one of the 
contacts is close to Gq. In this case, the G{Vg) function 
does not exhibit sharp peaks. Nevertheless, the question 
regarding the fluctuations of the distance between the 
conductance maxima is still valid. We study the origin 
of these fluctuations and estimate their variance. Fur- 
thermore, we extend the consideration to the region of 
intermediate values of the junction conductances. In this 
region, the peak-to-valley ratio of the G{Vg) dependence 
is large, but peaks are not sharp (their width exceeds 
the level spacing in the dot). Thus the developed theory 
yields understanding of the evolution of the peak spacing 
fluctuations with the strength of the dot-lead coupling, 
ranging from weak tunneling to weak reflection in the 
contacts between the dot and the leads. 

In this paper, we base our statistical description of 
the electron states in the dot on the random matrix the- 
ory. Despite the apparent failure of RMT in explaining 
the peak spacing fluctuations in a closed dot, this is an 
adequate approach for our purposes, because we concen- 
trate on the contribution to the peak spacing fluctuations 
coming from the non-zero conductance of the contacts. 
This contribution is determined by the randomness of the 
wave functions near the contacts. Fluctuating conduc- 
tance peak heights, which depend on the same pawdom 
quantity, apparently are well described by RMT .00 

In a partially open dot, the structure of the electron 
wave functions in the vicinity of the contacts affects both 
the amplitudes and positions of the conductance peaks. 
This results in a specific cross-correlation between these 
two characteristics of the random function GiVg). The 
correlation grows with the conductance of the contacts; 
we calculate the corresponding correlation function. 



II. CONDUCTANCE PEAK POSITIONS 

A. RMT Hamiltonian 

In our consideration, the quantum dot is described by 
the Hamiltonian of the constant interaction model: 



Ec 



k 



CiC 



k^k 



(1) 



In the framework of this model, the Coulomb repulsion 
of the electrons in the dot is accounted for by a relatively 
simple charging term [second term in (Q)], and is char- 
acterized by only one parameter, namely, the "charging 
energy" of the dot, Ec = e^/C, where C is the total ca- 
pacitance of the dot. The gate voltage Vg is represented 
by dimensionless parameter = CgVg/e, with Cg being 
the capacitance of the dot with respect to the gate. The 
(random) one-electron spectrum in the dot is represented 
by the first term in Eq. (^ , the energies Sk are measured 
from the bottom of the band. 



The weak dot-lead coupling can be accounted for by 
the tunneling Hamiltonian 



k,p 



^kp^\^p 



H.c. 



k,q 



^kq^k^^q 



H.C. 



(2) 



Here indices p and q denote the states in the left and 
right lead respectively. The dimensionless conductance 
of the dot-lead contacts is given by 



gL = (2^)2 (life, 



I ' A ' 



5« = (2^)^(|ife,|^)-^, 



(3) 



where (...) denotes the averaging over a statistical ensem- 
ble, A = {sn+i — £n) is the average level spacing in the 
dot, and Va (« = L, R) is the average density of states in 
the leads. A reflectionless dot-lead channel corresponds 
to g — 1; the conductance of such a channel equals Go 
(Go/2 for spinless fermions). 

The tunneling matrix elements are proportional to the 
values of the electron wave functions at the points of 
the dot-lead contacts. The randomness of the wave func- 
tions results in the fluctuations of the tunneling matrix 
elements in the Hamiltonian (||). The fluctuations of the 
wave functions in macroscopic leads can be neglected, 
and the tunneling matrix elements can be written in the 
form 



^kp 



ikq 




(4) 



In Sec. |T1-|T3 we consider the model case of spinless 
fermions, which is easier to follow. The proper extension 
onto t he ca se of real electrons with spins is performed 
in Sec. HD. Since at high temperatures the mesoscopic 



fluctuations are washed out, in this paper we concentrate 
on the low-temperature limit, T < A. 



B. Conductance peaks for spinless fermions 

Near a Coulomb blockade peak, where conductance 
through the dot is mediated by one resonant level, the 
transmission coefficient of the dot at the Fermi level is 
given by the Breit-Wigner formula 



Ll R. 



(5) 



where Fq, is the width of the resonant level with respect 
to electron tunneling to the a-th lead, and €{Af) is the 
deviation from the resonance. The center of the peak 
in the zero-bias conductance G(A/') for spinless fermions 
corresponds to such value of Af that e = 0. 
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For an almost isolated dot {tkp , t^q — 
near the n-th peak is given by 



0), function e(7V) 



(6) 



Here E^\j\f) is the ground state energy of an isolated 
dot with n particles [Hamiltonian (^] . The values of Af, 
at which the conductance has a peak are thus determined 
by the equation 



E. 



(0) 



1 (A/-i°^ 



(7) 



Therefore, for an almost isolated dot, when e is given by 
Eq. (pf), the spacing between the peaks is 



(0) 



1 



£n+l 



Ec 



(8) 



In this limit, the random component of the peak spac- 
ing has the same statistics as the level spacings A„ = 

£n+l — En- 

Upon the conductance increase, the electron transport 
through the dot is still dominated by the tunneling via 
the resonant level, as long as the effective level width re- 
mains small compared to the level spacing. The parame- 
ters entering in Eq. (^), however, acquire corrections due 
to the virtual processes involving other levels in the dot. 
The peak position is affected only by the corrections to 
e. At weak dot-lead coupling, the latter can be evaluated 
in the second order of the perturbation theory in tkp, tkq, 
so e is given by 



El^\Af)+6EUAf) 

-[£;i°2i(AA) + 5<_i(AA) 



(9) 




k=n+l 



ip + en-£k + Ec{n-Af-l/2) 
\tkp\'^0{-Cp) I 



ip-£n+£k + EciAf-n+3/2) 



(10) 



where {p —^ q} denotes a term identical to the first term 
in braces in Eq. (|l0|), but with all p's replaced with g's. 
Here SE'^ is the correction to the energy of the state with 
n electrons in the dot due to the admixture of all the 
states with n ± 1 electron, except for the lowest-energy 
state with n — 1 electrons. For spinless fermions, the mix- 
ing of the ground states with n and n — 1 particles in the 
dot determines the width of the n-th resonance (j|) but 
does not affect its position. The equation for 5E'^_i is 
similar to Eq. ® . 

These corrections to e shift the peak positions. 



(11) 



When the corrections 5E[-^{M) are small, the peak shifts 
9-^6 ^-Iso small and one can expand the right hand 
side of Eq. (|) in powers of 5J\f in the vicinity of Mn^ . 
Then it follows from Eqs. @) and (|) that e(7V) = when 



5NL 



SEL - SE' 



"SaT 



5E' - 5E' 



Er 



A/'=A/'<'" 



(12) 



Substituting the expressions for SE' i and SE'^ 
[Eq. (|l^)] into Eq. (|l2|), and using Eq. (|), wc find the 
shift of the conductance peak resulting from the dot-lead 
coupling: 



6K = 



1 A 

{27t)^E^ 



a=LM. k 

\M^o.)\\^^Ec + \e 



ek\ 



(13) 



The ensemble average of this shift equals zero. However, 
the mesoscopic fluctuations of the wave functions in the 
dot result in non-zero random shifts of the conductance 
peaks in each particular realization. Equation ( p^ yields 
the following expression for the fluctuating component of 
the spacing between the nth and in + l)st conductance 
peaks: 



5UL 



-■5K- 

A 

' E^ 



5ML 



E 



7a 



a = L,R 
/— n,n+l 



(14) 
(15) 



where 



(27r 



1 

In ^ 



(16) 



The term 5U'j^ in Eq. (|Tj) describes the contribution of 
the dot-lead coupling to the fluctuations of the spacing 
between the Coulomb blockade peaks. This contribution 
is statistically independent of the term Un^ , which ac- 
counts for the fluctuations in the energy spectrum of the 
dot and is given by Eq. (^). The logarithmic factor in 
Eq. (Ilj) is the result of summation over the lead states 
between the effective upper and lower cut-offs, Ec and 
A. The contributions of the wave functions of the states 
other than n and n -|- 1 do not contain the large factor 
ln(£^C'/A) and are neglected. 
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C. Spinless fermions: renormalization group 
approach 

Let us discuss the region of applicability of Eqs. (]l5|)- 
(p^). The right-hand side of Eq. (p^, obtained in the 
lowest non- vanishing order of the perturbation theory in 
tkp, tkq, is proportional to the product of the small pa- 
rameter ga and the large factor ln(i?c'/A). Examination 
of the higher orders of the perturbation theory shows that 
the leading term of the mth order is proportional to 



Ec 



Ec 



(17) 



Thus the results presented by Eqs. (15)-(jT^) are valid 
only if the condition 



/g^a-i — < 1 



(18) 



is satisfied. However, beyond the limits of the applica- 
bility of the finite order perturbation theory, which are 
given by Eq. ([T8|), one can exploit the leading logarithm 
approximation. It consists in summation of the leading 
terms [see Eq. (|l^] of aU orders of the perturbation the- 
ory. This can be done by means of the same renormal- 
ization group technique, which waainitially developed by 
Anderson for the Kondo problem.EJ Application of this 
method to the Coulomb blockade problenJl3 allows one 
to obtain results in the domain of parameters 



-1 < , 
,ln— <1 



(19) 



which is wider than the region of applicability (|18|) of the 
finite order perturbation theory. The technique used in 
Ref. as well as the original technique by Anderson, 
does not account for the mesoscopic fluctuations of the 
one-electron wave functions. For our purposes, we need 
to extend the renormalization group method to account 
for the randomness of the electron states in the dot. Be- 
low we present the proper generalization valid for a dot 
obeying the RMT statistics. 

The Anderson's scaling procedure consists in trans- 
formation of the initial Hamiltonian Hq + i/t to a new 
one, with smaller bandwidth and with the matrix ele- 
ments renormalized to compensate for this band reduc- 
tion. The transformation operates within the invariant 
space of Hamiltonians having generic form 



E 



zU) 



rilV2 Vl V2 



mm 



P. 



(20) 



where the indices rji , 772 run through the one-particle 
states in the dot and both leads, and Pj is the opera- 
tor of projection to the Hilbert subspace of states with 
j particles in the dot. Near the n-th peak, the renor- 
malization of the Hamiltonian occurs mostly due to the 
virtual transitions between the states with n and ?i — 1 



particles in the dot, and it is sufficient to limit our con- 
sideration to these two sets of states. Then the summa- 
tion in Eq. (|2(j) must be done over j = n — l,n only. 



and T 



(") 

kp 



t(") 

' kq ' 



-p(n- 
' qk 



= because the corre- 



sponding transitions involve states outside this subspace 
of states. 

The transformation starts from Hamiltonian ([|)-(||) 
operating within the electron band of width Dq ^ Ec- 
In the course of renormalization, the band width D is 
reduced in steps starting from D = Dq. As we men- 
tioned above, the matrix elements must be adjusted to 
compensate for the band reduction: after each step, the 
transition amplitudes calculated with the renormalized 
Hamiltonian H(Z3 — dD) must coincide with those de- 
rived from the Hamiltonian H(D) of the previous step. 
In the leading logarithmic approximation, which we em- 
ploy, these amplitudes should be calculated in the second 
order of the perturbation theory. Then the scaling law, 
which dictates how a matrix element must be changed 
when the bandwidth is reduced from D to D — dD, in 
general reads 



dT^^^ = 



E 



' Vim ' mv2 



D-dD<\i„^\<D 



V3\ 



(21) 



where £[r]2 — > rj^] is the energy difference associated with 
the (virtual) transition 772 — >■ ?73, and js is the number of 
electrons in the dot after this transition. Here we present 



(") . 



explicitly only the equation for T^^^ 

(n) \ ^ "Tfc^pTpfc^ 



dJ 



-D<^p<-D+dD 



j(n) j(n) 



D 



E 

-D<^g<-D+dD 



Tfc, „T, 



kiq ' qk2 



D 



- E ■ 

D-dD<^k-i<D 



D 



- E ■ 

-D<ik^<-DJrdD 



-p(n) -p(n) 



D 



(22) 



The equations for the other parameters of the Hamilto- 
nian have similar form. The renormalization procedure 
also generates the constant term Eq: 



E 



' mm ' mm 



-D<ini<-D + dD 
D - dD < C„ < D 



2D 



E 

-D<i^<-D+dD 



mn 



(23) 



The initial Hamihonian H(D = Ec) = Ho+Ht [Eqs. @)- 
(^] in terms of Eq. (20) with j = n — l,n corresponds 
to 



T 



kp 



0, 

T* 
' pk 



tkp 5 ~^ kq 



T 



qk 



^kq •) 



(24) 



tu) - T^-J') - n t(^') - 



It is random because of the fiuctuations of one-electron 
wave functions [see Eq. (0)] and level spacings. The 
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Hamiltonian H{D < Ec), which is produced by the scal- 
ing transformation from Hq + Ht, is therefore also ran- 
dom. 

As it can be seen from Eq. (p^), the randomness of the 
initial Hamiltonian Hq + Ht is inherited by the renor- 
malizcd Hamiltonian H(Z?) in two ways. First, there are 
random contributions from the states fci and k2. Second, 
there are also random contributions from the states /ca, 
being eliminated by the band reduction. The final renor- 
malized matrix element Tfe^fc^ "^iH still be proportional 
to the (random) wave functions V'^j and V'fej- C)n the 
contrary, the random contributions from the eliminated 
states fca are added together and their (independent) fluc- 
tuations are averaged out. This consideration leads one 
to the conclusion that the solutions of the RG equations 
( pl| ) can be cast in the form: 



±A 



E 

a=L,R 



V'fci(l-a)V'fc,(ra) 



77=^ /"il^i : 



T—\r{D) 



X^iD) , (25a) 

(25b) 

(25c) 

(25d) 
(25e) 



Here the upper sign in the right-hand-sides corresponds 
to J = n — 1 and the lower sig n corresponds to j — n. 
From Eqs. (|2l|), (^J), and (^, we derive the scahng 
equations for the non- fluctuating parameters and fii, 
which determine the behavior of the Hamiltonian in the 
course of the scaling procedure: 



dlnD " d\nD 

with the initial conditions 



4X„{D)^la,{D) (26) 



The transformation can proceed until the band width 
is reduced to the value of the order of the mean level 
spacing A. After the renormalization, the ground state 
energy can be found with the help of Hamiltonian H (D ~ 
A) as the sum of Eq and the correction due to the dot- 
lead tunneling in the reduced band. As we explained 
above, the latter term is the main source of fluctuations 
of the ground state energy. It can be calculated in the 
second order of the perturbation theory in T^^^^, simi- 
larly to Eq. ([lO|). This calculation finally yields Eq. ( |T5| ) 
with 



7'=-- 



/^tan 



1 Ec 

TT A 



(29) 



If the argument of the tangent function is much less than 
unity, Eq. ( p9| ) is reduced to Eq. (|l^), as expected. 



D. Electrons with spin 



In this section we will refer to the same Hamiltonian 
(|l|)-(H), with summation over the spin index added. The 
conductance peaks are numbered as follows: the (2to— 1)- 
th conductance peak corresponds to the resonant tunnel- 
ing of an electron through the empty m-th level, and the 
(2TO)-th peak corresponds to resonant electron tunneling 
via the m-th level already occupied by one electron. 

In principle, the spin degeneracy can lead to formation 
of the many-body state, which enhances the conductance 
in the valley between the {2m~ l)-th and (2m)-th peaks, 
so Eq. (^) is no longer valid. This phenomenon is com- 
monly referred to as the Kondo effect. However if the 
temperature is larger than the widths of the dot levels, 
which is the case under consideration, the formation of 
this many-body state is suppressed, and we can employ 
the rate equation formalism to evaluate the conductance 
of the dot. This approach^ields the following expression 
for the dot conduct ance:t2l 



X^{D^Ec)^Q. fic.{D = Ec) ^ —Vg^ . (27) 

They coincide with those derived in Ref. |l^ and have 
solutions 



47r 



^J■a{D) = -^Vda i COS 



1 Ec 
n U 

1 ^, Ec 

TT D 



(28) 



The energy shift Eq also accumulates in the course of 
renormalization. Its random part, which stems from the 
fluctuations of the wave functions of the eliminated states 
A;3, averages out due to the summation over a large num- 
ber of these intermediate states, similarly to the corre- 
sponding contribution to the matrix elements Tp^p^ and 



G{Af) = Go 



9L9L dfF[e{Af)]/de 
gL + gRl + fF[±eiM)] 



(30) 



where /f(0 = l/[cxp(C/^) + 1] is the Fermi distribution 
function. The sign in the denominator should be taken 
for even peaks and "— " for odd peaks. One can 
see that the peaks in the conductance G{J\f) are shifted 
from the points of charge degeneracy, which correspond 
to e(A^) = 0, by ±(ln2/2)T, which results in the correc- 
tion to the peak spacings: 



SU 



(T) 
2m- 1 



(In 2) 



T 

E^ 



5Ul 



(T) 



2m 



-(In 2) 



Er 



(31) 



For an almost isolated dot, the distances between the 
charge degeneracy points are given by 



5 



(0) 

(0) 
2m 



u 



1 



(32a) 
(32b) 



To obtain the expressions for the conductance peak 
spacings for a partiaUy open dot, one has to add up U'^'^\ 
5U'^^\ and the contribution to the peak spacing due to 
the finite dot-lead tunneUng. This contribution has two 
components. 

The first com pone nt, 51A'^, is analogous to the one cal- 
culated in Sec. II B, it comes from the transitions in- 
volving all levels in the dot except for the resonant one. 
The spacing U2m is between the conductance peaks corre- 
sponding to filling of two consecutive orbital levels. The 
result for spinless fermions is given by Eq. (15). For elec- 



trons with spin the corresponding correction acquires an 
additional factor of 2: 



J A 



E 



7a 



1 - 



a = L,R 
/— m,m+l 



IV'/(r.)P 



(33a) 



The correction of this type contribute equally to the po- 
sitions A/'2m-i, J^2m Or the peaks, corresponding to the 
same orbital level. As the result 



0. 



(33b) 



The other component, dlA'^^ is due to the virtual transi- 
tions from/to the resonant level. Let us consider, for def- 
initeness, the 2m-th conductance peak, whose position is 
determined by the energies of the ground states |2m — 1) 
and 1 2m), with 2m— 1 and 2m electrons in the dot respec- 
tively. For spinless fermions, virtual transitions involving 
the m-th level do not lead to the shift of the peak po- 
sition, because of the electron-hole symmetry. The spin 
degeneracy, however, removes this symmetry and as the 
result the shifts of the-ffround state energies 5E2m and 
(5i?2in-i are not equahliil The state |2m) acquires an ad- 
mixture of the excited states with 2m — 1 electrons in the 
dot, and the corresponding correction is given by 



. P 



I ^mp 



l^mqP^('Cg) 



(34a) 



where the prefactor 2 is due to the spin degeneracy. The 
state 1 2m — 1) intermixes with the excited states having 
2m or 2m — 2 electrons in the dot. At the point of the 
2m-th conductance peak, which we consider, the states 
|2m — 1) and |2m) have equal energy, but the energy 
of the state |2m — 2) is higher, because of the charg- 
ing. Therefore the coupling of the states with 2m — 1 
and 2m — 2 electrons in the dot is suppressed and the 
correction 5ii'2m-i is smaller than SE!^^: 



E 



E 



I ^mp 

^p + 2Ec 



(34b) 



Substitution of Eqs. (§4) into Eq. (ha) yields 



|V'm(ra)|^ 



with 



1 1 



(35) 



(36) 



The shift of an odd peak can be calculated analogously 
and is given by 



(37) 



In contrast to the correction from the non-resonant lev- 
els [Eq. (p^)], the correction 5Af^ has a non-zero ensemble 
average, 



a=L,R ^ 

Equations (H), @i yield 



(38) 



'5^2™ = Tr; 2. 7a- 



Ec ^ m 



(39a) 



a = L.R 
Z— m,m-t-l 



J'/'m(ra^_ (39b) 



m 



The full expression for the peak spacings in the case of 
electrons with spin is 



Ur, = + m^P + 5U' + SW' 



(40) 



where SuP and U^"^ are given by Eqs. (|l]) and 1^) 
respectively. 

Equations (|35|)-(|3^), derived with the help of the sec- 
ond order perturbation theory in tkp, tkq, are valid only 
at relatively small values of the conductances ga of the 
contacts (see Sec. IIC ). When the perturbation theory 
fails, the proper expression for 7^' can be obtained with 
the help of the renormalization group method, described 
in Sec. [IC. If the temperature is of the order of the 



mean level spacing, then the RG transformation must 
be stopped when the bandwidth D is reduced to A; the 
renormalized 7^' equals 7^, which is given by Eq. (p^). 

If T < A, then we have to continue the renormaliza- 
tion procedure to the bandwidths D < A. There is only 
one level in the dot within such a band, so Eqs. (^), ([25| ) 
yield: 



d\nD 



= Ma(^) 



dfJ-ajP) 

dlnD 



= 2fi^{D)X^{A) . (41) 
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The sohitions of these equations determine the renor- 
mahzation of the tunneh ng m a trix elements Tkik2i T^p, 
and Tkq given by Eqs. (25a)-(25c). As for Tp^p^ and 



Tqjqj, th ey ar e not renormalized anymore and are given 



by Eqs. (25d)-(25e) with D — A. The renormalization 
procedure should be carried out until D ^ T. Finally we 
extend our result ( ^sf ) to a wider region of parameters, 
defined by Eq. ([l9|), with parameter 7^ given by 



7" = 7' 

la I a 



2A„(A) 



2A<:,(A) 



- 1 



(42) 



Here Aq(A) and /in(A) are determined by Eqs. (Pq). The 
first term in Eq. ( [42| ) is produced in the course of renor- 
malization while the band is reduced to A; the second 
term is produced by the band reduction from A to T. 
Like in the spinless case, Eq. (|4|) reduces to Eq. ( ^ ) at 
small values of the dot-lead conductance, i.e. when the 
condition (nm is satisfied. 



III. STATISTICAL PROPERTIES OF THE PEAK 
SPACING 

A. Peak spacing fluctuations 

Now we discuss how the finite tunneling strength af- 
fects the distribution function of the peak spacings, which 
is given by 



P„{U) = / dW(°)p(°)(MW) Y[ dM(ra)P({V'}(r„)) 

a=L,R 

X 5{U-Un) ■ (43) 



Here Un is determined by Eq. (|4C|). It follows from 
Eq. (H) that the distribution function oi U^"\ P'"', is 
a (5-function for n — 2m — 1, and is given by the Wigner- 
Dyson distribution function for n = 2m. The set UjA 
denotes the wave functions entering Eqs. (|3^) and ( |39[ ) 
for 5U'^ and 5U'^. The values '(/'m(r) are distributed ac- 
cording to the Porter-Thomas statistics. 

At small junction conductance, corrections SW^, SU^ 
are negligible. In this limit, the distribution of the peak 
spacings is manifestly bimodal. Odd values of n corre- 
spond to the spacing between two peaks originating from 
filling of the same orbital level in the dot; the correspond- 
ing distribution function is infinitely sharp, and is cen- 
tered a.t U = 1-|-Tln2/ Ec ■ Peak spacings at even n 
correspond to filling of the consecutive orbital levels, and 
therefore the distribution function for them follows the 
Wigner-Dyson statistics. Note, that at a finite temper- 
ature the centers of the "even" and "odd" distributions 
are shifted towards each other, so these two distributions 
overlap. This shift, which is the result of interaction be- 
tween the electrons tunneling via a spin-degenerate level, 
was not taken into account in the analysis of experimen- 
tal data.Ll 



At finite contact conductance, mesoscopic fluctuations 
of the wave functions in the dot yield a finite width of the 
"odd" distribution and broaden the "even" distribution. 
These effects stem from a finite value of 6U^^ and SU^. 
The resulting widths are characterized by 



{{SU 



2m- 1) 



Eh 



(44b) 



Here /3 = 1 and 2 for orthogonal and unitary ensembles 
respectively. For smaller values of junction conductances 
ga [such that condition ( p^ ) is satisfied], ^'^ and 7^' are 
proportional to ga and are given by Eqs. ( p6| ) and ( |3^ ) 
respectively. Larger conductances cause renormalization 
of 7^ and 7"; the resulting expressions for them are given 
by Eqs. @) and (||). 



One can see from Eq. ( 441 ) that there are two indepen- 
dent additive contributions to the fluctuations of the peak 
spacing. The first one is due to the fluctuations of the 
wave functions and is proportional to the junction con- 
ductance. The second one comes from the fiuctuations of 
the electron spectrum. If the latter is described by RMT, 
then {SAD = [(4/7r) - 1] A^ and {SA^ = [(Stt/S) - 1]A^ 
for orthogonal and unitacy ensembles respectively. The 
experimentally measuredlH {SA^ apparently are larger 
than the values predicted by RMT. Regardless this dis- 
crepancy, the contribution coming from the randomness 
of the wave functions can be singled out from experimen- 
tal data by comparing the values of {SU^ measured at 
various junction conductances. 

The average spacing between the peaks is also affected 
by finite conductance: 



(Z^2™-i) = l + (ln2)-|^-2 ^'^4-^ 

^ a=L,R ^ 



(45a) 



(W2») = l + ^-(ln2)^ + 2 E 7:A. (45b) 



One can also consider the distribution function aver- 
aged over even and odd peak spacings. A sample plot 
of this function is shown in Fig|-Jl|. It is instructive, fol- 
lowing the experimental paper ,0 to compare this func- 
tion to the distribution function P^^^'ipl) of the spacings 
between the charge degeneracy points in an almost iso- 
lated dot. The "bimodal" distribution P'^^^iJA) consists 
of the Wigner-Dyson distribution ( "even" spacings) and 
a (5-peak ("odd" spacings). The key differences of P(U) 
from P(°)(iY) are as follows: (i) the 5-peak corresponding 
to odd spacings is broadened (ii) the centers of the dis- 
tribution maxima corresponding to even and odd peak 
spacings are shifted by ±[(ln2)r/£;c - 2(^^ ia)^/Ec] 
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(iii) unlike P^°-' {U)^ the distribution P(ZY) has a tail which 
extends into the region U < \. 




12 3 

(W„-1)/(A/Ec) 

FIG. 1. The cumulative distribution function of the "even" 
and "odd" spacings between the Coulomb blockade peaks, 
W„/(A/£;c), given by Eqs. (g), (||) for 7^ = 0.075, 7^.' = 0.1, 
T = 0.5 A (solid line). The dotted line shows the "bimodal" 
distribution P'°' consisting of the Wigner-Dyson distribution 
( "even" spacings) and a (5-peak ( "odd" spacings) . The shown 
distributions correspond to the Gaussian unitary ensemble. 

One can see from Eq. that the fluctuations of wave 
functions and of the level spacings contribute equally to 
the width of PiU) at 7^ ~ 7^' ~ 0.15. The dimensionless 
parameters 7^ and 7^' depend not only on the junction 
conductances, but also on the ratios Eq / ^ and Ec/T, 
see Eqs. ^ and @. Using the values Ec = SOOueF, 
A — 7/iey , and T ^ A from the experiment of Ref. 0, we 
obtain that the values 7^ ~ 0.05 are reached when the 
junction conductances equal to O.SGq. Thus a consid- 
erable contribution of wave function fluctuations to the 
randomness of the interpeak spacings can be achieved 
when the conductance peaks are still well-defined. 



X n dW(ra)P(W(r.)) 

a=LM 

x<5(zi-Z^„) ^(g^'"^'^) -Gi'"'^'')) . (47) 
Here the spacing of the n-th peak from its neighbors is 

Un = 2 ' 

and Un is given by Eq. (^). The distribution 
P{U,G^^'''''^) is skewed, i.e. {SU SG'-'^''''^) ^ 0, because 
of the dependence of U on {-0}. The dimensionless cross- 
correlation parameter grows with junction conductance, 



Kg-u 



(SG 



(max) 



{Gir^'M^/Ec) 

- E ^(f^)(27;+7^:) 

a' ^ a 



(48) 



Function J-{x) > is given by 



2Vi 



l + Vx 

for the orthogonal ensemble and by 
2x^ + 3x^ - 6x + l- 



T{x) 



(a;-l)(x2-l-2a;lnx) 



for the unitary ensemble. In the case ql — Qr, the argu- 
ment of the function T in Eq. ( p8| ) equals unity; = 1 
and 3/2 for the orthogonal and unitary ensemble respec- 
tively. 



B. Correlations between the heights of the peaks 
and their spacings 

We have shown that for a partially open dot, the con- 
ductance peak positions depend on the values of the one- 
electron wave functions at the points of the dot-lead con- 
tacts. On the other hand, these values determine also 
the conductance peak amplitudes: 

^2ra-l-^2^ °^ 5^ |0„ (r^) 1 2 + |^„^ (r^ ) 1 2 " ^^Oj 

Therefore the spacing between two peaks must be cor- 
related with the heights of these peaks: the higher the 
peaks are, the less (on average) is the spacing between 
them. 

The joint distribution function of the conductance 
peak height G^"^'^^^ and the spacing of this peak from 
its neighbors U is given by 



C. Fluctuations of spacing between distant peaks 

In the previous sections we considered only the fluctu- 
ations of spacing between the neighboring peaks. Now 
we will discuss the statistical properties of the spacing 
between the n-th and (n -I- Z)-th peaks, 

Un{l) = Nn+l - Nn , (49) 

with ^ > 1. Here we must mention that in real dots 
adding-electrons apparently scrambles the onejelec.tron 
states.EJ However, according to the experiment ,EjI1j the 
changes in the dot's one-electron spectrum and wave 
functions are small when the number of added electrons 
is small. Under these conditions, the same realization of 
the random parameters of the Hamiltonian (l)-(2) de- 
termines the positions of ?ith and [n + ^)th peaks. The 



8 



average spacing between the peaks can be derived from 
Eqs. (|3l|), and (|38|), and is given by 



{Unil)) =l[l + 



A 

2Ec 



A 



1 ln2 T 

2 ^ ~A 



[(-1)"+' - (-1) 

E 7;: 



a=L,R 



(50) 



Like in the case of neighboring peaks, which we consid- 
ered in Sec. |l[ the fluctuations of the peak spacings have 
the independent contributions SW , 6U" and U'^^^ com- 
ing from the randomness of the wave functions and level 
spacings respectively. Using Eqs. (p^), and (^), we 
obtain for / < £'c/A 



-V 



[4(Z-l)(7if + , (51) 



a=LM 



where depend only on whether the numbers n and 
n -|- I are even or odd, and do not grow with I: 



Ba — 



ill? 

(27; 

ill? 



even I , 
, odd Z > 3, odd n 
I = 1, odd n , 
odd L even n . 



Thus the contribution due to the fluctuations of the wave 
functions, 6U^{l)+6U^{l), is proportional to the junction 
conductance, and its variance grows linearly with the dis- 
tance I between the peaks. This behavior is in qualita- 
tive agreement with the experiments.E3 The fluctuations 
of the other contribution, SUn\l), do not depend on the 
junction conductance, and are characterized by 



' ^fi + i-ir+'-i-iy 



(52) 



Function F(k) determines the variance of the spacing be- 
tween two distant one-electron energy levels, and is given 
by the RMT; F{k) In/c at A: > 1. One can see that 
for sufficiently distant peaks and large junction conduc- 
tances, 5U'{1) + 5U"{1) dominates over 5U^^\l), and the 
main term in the fluctuations of the spacing comes from 
the fluctuations of the wave functions at the points of 
dot-lead contacts. 

At large the correlations between the sets of the lev- 
els that determine the shifts of the n-th and {n + 
peaks become negligible. If the Hamiltonian of the dot 
is not aS^ected by adding of new electrons, it occurs at 
I ^ Eel ^1 when these two sets of levels, having char- 
acteristic width Ec each, no more overlap. In real dots, 
it can occur even at smaller Z, because by adding elec- 
trons one modifles the effective potential which conflnes 



electrons to the dot. Eventually the modifications result 
in a loss of correlations between the one-electron states 
in the original and new potentials. Then the positions 
of distant peaks are defined by two Hamiltonians, both 
of the generic form ([l])-(H) with the same Ec^ but with 
different and statistically independent (random) t^j,, t^q 
and £fc. The resulting variance of the peak spacing at 
large / is twice the variance of a peak position [Eq. (|l2|)] 
and is given by 



1 



67r2/3 Ec 



E 

a=L,B. 



9a 



(53) 



The contribution Shl^{l) does not grow with I and is ne- 
glected here. 

Thus the fluctuations of the interpeak spacing due to 
randomness of the wave functions saturate at large dis- 
tances between the peaks. The corresponding fluctuation 
amplitude exceeds that for the neighboring peaks by a 
parametrically large factor ^/EcJA./[Iii{Ec / A)]"^ ■ This 
saturation of the peak spacing fluctuations with the in- 
crease of I is in agreement with the experimental results 
(see Fig. 2 in Ref. |l|). 

As for the contribution of the level spacing fluctua- 
tions, Un\l), its statistical properties for large I de- 
pend crucially on the nature of the scrambling of the 
one-electron states, and are model-dependent; the dis- 
cussion of this issue is beyond the scope of the present 
paper. The experimental verification of Eq. (|5^), how- 
ever, is possible without knowledge of {[5lAn\l)Y)- Since 
5lAn\l) is independent on the junction conductance, the 
difference between the values of {[5Un{l)Y) measured at 
various junction conductances can be directly compared 
to Eq. (||). 

Equations (|5l| ) and (^3|) describe the behavior of 
{5Uf^{l)) in the limiting cases of small and large I re- 
spectively. The crossover between these two limits can 
be described analytically in the absence of scrambling of 
the electron wave functions by adding new electrons to 
the dot. Equations (fsl), (pll) and (p5|) yield 



32 



a \Ecj 



E 

a=L,R 



(54) 

where Xa{D) is given by Eq. (p8|). At small /, one can 
replace Xa{D) in Eq. ( [5^ ) with Aq(A) — 7^, which yields 
the Z-dependent part of Eq. (|5l]). For / ^ Ec/A, one 
can replace Xa{D) with ga:/{2TFp within the most of the 
integration interval, so the right-hand-side of Eq. ( [5^ ) 
matches Eq. (pSh. 



IV. 



CROSSOVER BETWEEN WEAK AND 
STRONG TUNNELING 



In the previous sections we considered the case of small 
enough conductances of both dot-lead contacts [the con- 
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dition (|T^) is satisfied]. In this weak tunneling regime, 
the n-th conductance peak corresponds to the resonant 
tunnehng via the single, n-th level. The width of the 
peak is given by 



En 



(55) 



where Tna is the partial width of the n-th level with re- 
spect to the electron tunneling to the a-th lead. If the 
coupling of the dot to the leads is weak enough so the 
condition ( [l8| ) is satisfied, the partial level widths are 
given by the conventional expression 



|Vn(ra)P A 



2tt 



(56) 



At larger values of ^/ g^ In jEg/A), the level widths are 
renormalized (see Sec. |ll C| ): 



r„„ = [27r/x„(A)] 



,| V>„(r^)P A 
(l^nP) 27: 



(57) 



It implies the renormalization of the junction conduc- 
tance: 



(58) 



where ^ia{D) given by Eq. (psj). 

The position of a conductance peak is mainly deter- 
mined by the position of the corresponding charge de- 
generacy point. However, there is another random con- 
tribution to the shift of the peak position, which stems 
from the randomness of the phases of the transmission 
amplitudes via different one-electron levels in the dot. 
Even when the dot is in the resonance, and the dominat- 
ing contribution to the tunneling amplitude comes from 
a single resonant-level, there is also a contribution from 
the cotunnelingE3 through the other levels, which yields 
the following estimate for the tunneling amplitude: 



[{T,,L+rnR)/2]^ + {Ec6Afr 



vT/cLTfe 



kR 



£k — £r 



EcSM 



(59) 



where 77^ is the random phase of the amplitude of tun- 
neling through the k-th level in the dot, and SAf is the 
deviation from the closest charge degeneracy point. The 



position of the maximum of |yl„((SA/')| given by Eq. (59), 
SAf = SAfn, may not coincide with the point SAf = 0, i.e. 
the cotunneling contribution adds to the random shift of 
a peak. Equation (^ , together with Eq. ( |57| ) , yields the 
following estimate for the corresponding contribution to 
the peak spacing fluctuations: 

^^(cotunn) ^ ^ [^2 ^ ^2^(^)]3 ^^^^ 



This contribution is much smaller than the one com- 
ing from the fluctuations of the charge degeneracy points 
[Eqs. (H), (HI)], provided condition (|l9|) is satisfied. This 
allowed us to neglect the contribution ( |60| ) yet so far. 

Now we would like to extend our consideration beyond 
the limits set by Eq. ( p^ . To demonstrate the new phys- 
ical features emerging at stronger coupling of the dot to 
the leads, it is sufficient to limit ourselves to the case 
when the conductance of one channel (say, the right one) 
is small, and satisfies Eq. (p^), as the other one (5^) 
varies between and 1. The consideration of the whole 
2-dimensional parameter space < < 1, < g/j < 1 
would not add anything new to our understanding of the 
nature of the fluctuations of the conductance peak posi- 
tions and will be omitted. 

To begin with, we consider the conductance of the dot 
in the conditions opposite to the limit of weak coupling. 



5fl < 1 , 1 - 5L < 1 • 



(61) 



This case was studied in details in Ref. |l^ and we first re- 
call some findings of this work. Since one of the contacts 
is almost open, the charge of the dot is not quantized 
anymore. The Coulomb blockade is almost lifted and the 
notion of the charge degeneracy points is useless. In- 
stead of sharp peaks, the conductance as a function of 
gate voltage Af exhibits oscillations with the period of 
unity. The oscillations occur around the value 



(G) = GagR 



_A_ 



(62) 



the characteristic amplitude of the oscillations is also 
GagR{A/ Ec). The phase of oscillations is a random 
quantity, varying with Af slowly (at scales larger that 
1). The electron transport through the dot is mediated 
not by a single resonant level, as in the case of weak cou- 
pling, but by a set of levels, whose energies lie in the strip 
of the width Ec near the Fermi energy. As a result, the 
value of conductance at given gate voltage Af is deter- 
mined by the properties of all levels in this set. Since 
there are no well-pronounced peaks, it is more natural 
to describe the statistical properties of conductance not 
in terms of the heights and positions of the conductance 
maxima, but rather in terms of the correlation functions 
(G(M)G(AA2)). 

However, we still would like to study the statisti- 
cal behavior|-e|f the conductance maxima, following the 
experimentsjlj and in order to make a connection to the 
case of weak conductance. The quantitative considera- 
tion is hampered by the fact that, in the absence of sharp 
peaks, the exact positions of the conductance maxima 
are not associated with special points of any physical 
quantity. Therefore we will limit ourselves to a qual- 
itative consideration of the fluctuations. It was men- 
tioned above, that the Coulomb blockade is lifted when 
f/L « 1. Under these conditions, the conductance max- 
ima occur at the points where the amplitudes of tunneling 
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via En/ A levels in the dot add up most favorably [Cf. 
Eqs. (js^), (|60|)]. The phase of transmission via each level, 
rjk, depends on the gate voltage, so the function G{J\f) is 
not constant. It is (almost) periodic, because the func- 
tions rjk{N) have period of unity; G{JV) is not exactly 
periodic because when the gate voltage Af is shifted by 
unity, one level is replaced in the set which mediates con- 
ductance. Therefore, the conductance maxima are not 
equally spaced, but their positions fluctuate with respect 
to each other. When Af is varied by Ec/ A, all the "con- 
ducting" levels are replaced, so G{J\f) is not correlated 
with G{J\f + Ec/ A). In particular, the positions of the 
maxima separated by Ec / A other maxima fluctuate in- 
dependently within range 6Afn ^ 1 each. Therefore, the 
characteristic amplitude of fluctuations of the spacing be- 
tween nearest maxima equals 



(63) 



for the conditions defined in Eq. (|6l|). 

Now we can combine the knowledge acquired in the 
consideration of the case of weak [Eq. (p^)] and strong 
[Eq. (|6l|)] tunneling through the left contact, to deal with 
the case of the intermediate conductance 




HEc/A) 



< 5L < 1 



(64) 



In considerations of Sec. ^ we used the renormalization 
group transformation of the Hamiltonian only in the very 
vicinity of the conductance peak. In principle, this trans- 
formation can be performed at any value of 6J\f, where 
SM is the distance to the closest charge degeneracy point 
of Hamiltonian (|l|) . If 



\SAf\ > exp(-7r/V5L) 



(65) 



then the renormalization must be stopped when the re- 
duced band width D reaches Ec\SJ\f\. Indeed, if D < 
Ec\SJ\f\, then one of the two dot's charge states is left 
beyond the band limits; only one charge state is then 
possible. One can see that the renormalization defined 
by Eq. (|2^) does not bring the system out of the condi- 
tions of weak couphng, because gL{D = Ec\SJ\f\) is still 
much less then unity. The conductance of the dot can 
then be estimated from H{D — Ec\SM\) with the help of 
the formula for elastic cotunncling. Under the condition 
(|65|), the conductance is 



G^GogngLiD)- , with D = Ec\SAf\ . (66) 

The renormalization of the conductance of the right con- 
tact, gR, can be neglected, since it satisfies Eq. (p^. 
If, on the contrary. 



then the renormalization is stopped when D reaches 
i?c exp(— tt/^/^l). At this point, the renormalized con- 
ductance of the left contact, given by Eq. (^81) , becomes 
of the order of unity. Then the scaling law ( pl[ ) , based on 
the perturbation theory in the small parameter gL{D), is 
not applicable. The nature of conductance through the 
dot with such parameters (one contact is almost open, 
the other one is almost closed) was discussed earlier in 
this section. The characteristic conductance for the con- 
ditions of Eq. (|65| ) is given by the formula similar to 
Eq. (|6^). In this formula, Ec, which in Eq. ( |6^ ) plays the 
role of the width of the band of electrons affected by the 
charging, must be replaced with D — Ec exp(— Tr/y^g^): 



A / vr 
G Gog R— exp 
Ec 



9l 



(68) 



This value is parametrically larger than the one given by 
Eq. (^. Thus the characteristic width of a conductance 
peak in the conditions of Eq. (pif) is 



>V„ ~ exp 



(69) 



The conductance maxima are not tied to the charge 
degeneracy points. Instead, their positions determined 
by the interference of the amplitudes of transmission 
through {Ec /A) exp{—n / y/gZ) levels, that belong to the 
reduced band of width Ec exp(— Tr/^/gi). Therefore, the 
conductance may have a maximum, roughly speaking, 
at any point lying in the interval where the renormal- 
ized conductance of the left contact is close to unity, i.e., 
within the region defined by \SJ\f\ < exp{—n/,JgZ)- The 
positions of the maxima corresponding to two neighbor- 
ing charge degeneracy points are correlated, because they 
are defined by the same set of levels except for one. This 
situation is similar to that in the case when one of the 
contacts initially had conductance close to unity. The dif- 
ference is that now the set of the levels which determines 
the position of a maximum is smaller. Employing the 
same ideas we used for the dot with one contact initially 
almost open, we obtain that the characteristic amplitude 
of fluctuations of the maxima spacings is 



(70) 




'9L 



\6M\ < exp(-7r/v^) 



(67) 



for the conditions defined by Eq. (|64|). This expression 
matches Eqs. (||) and (H) at gL ^ [tt / \n{Ec / A)]^ and 
gL ~ 1, respectively. 

Equation (|7^) gives an estimate for the fluctuations 
of spacing between two neighboring maxima. Now we 
extend this estimate to distant maxima. As we dis- 
cussed before in this Section, the random positions of 
the conductance maxima are determined by the sets of 
^ [Ec/ A) exp(— tt/y^) levels. For the n-th and {n + l)- 
th maxima, the levels in these two sets are the same ex- 
cept for I of them. These / levels in each of the two 
sets produce uncorrelated contributions to the positions 
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of the two maxima. This contribution grows with /, since 
the contribution of each particular level is independent 
of the others: 



SUnil) 




exp 



(71) 



for I < 



Er 



■ exp 



'9l 



At I > {Eel ^) exp{—n/y/^) there are no common lev- 
els in the two sets, and the fluctuations of the positions 
of the n-th and (n + l)-th peaks are independent. The 
variance of spacing between them saturates at large I: 



(72) 



SUnil) ~ exp - 



/9L 



for 



i > — exp 



The levels that are eliminated by the reduction of the 
band down to Ec exp^^n / ^^gZ) in the course of renor- 
malization also make a contribution to SUn{l)- However, 
it is important only in the narrow region of crossover 
between the regimes of weak [Eqs. (p^, (p3|)] and inter- 
mediate [Eqs. ([TI]), ([72|)] conductance, where the number 
of levels in the reduced band is of the order of unity. Oth- 
erwise this contribution is small, as compared to the one 
given by Eq. ([7l|), and can be neglected. 



V. CONCLUSION 



Coulomb blockaded quantum dot. We found that the 
fluctuations of the electron wave functions in the dot 
contribute to the fluctuations of the peak spacings; this 
contribution grows with the conductance of the contacts, 
see Table |. 

For relatively weak coupling of the dot to the leads, 
9L.R < [7''/ln(i?c/A)]^, the peak conductance is medi- 
ated by a single electron level, and the peaks are sharp. 
The variance of the peak spacing fluctuations for the 
neighboring peaks is given by Eqs. (^); it reaches 
when gL,R ~ [tt / \n{Ec / A)]^ . 

At larger values of gL,R, the conductance peaks are 
broad, and many levels are involved in the transport 
through the dot. Thus the properties of conductance 
peaks (more precisely, maxima) are determined by the 
wave functions of many electron states. The variance of 
the spacing between the maxima continues to grow with 
the conductance of contacts, and is given by Eq. ([70|). 

The fluctuations of the spacing between two conduc- 
tance peaks grow with the distance between them, see 
Eqs. (pl|), ([7l|). At large distances between two peaks, 
the variance of the spacing between them saturates at 
the value given by Eq. (H) (||). 

The dependence of the peak position on the one- 
electron wave functions in the dot leads to the corre- 
lations between the fluctuations in the peaks' spacings 
and heights: the spacings between the higher peaks are, 
on average, smaller. The parameter characterizing this 
correlation is given by Eq. (Esl). 



We have studied the statistical properties of the spac- 
ings between the conductance peaks in a partially open 



TABLE I. The characteristic values of the peak width and the amplitude of fluctuations of the spacing between the neigh- 
boring peaks in units of Ec for gu <^ 1, gn < ql < 1, as given by Eqs. (^), (|57|), (p9[), ([7o|). The numeric factors are 
omitted. 



conductance of the contacts 



peak width 



gR< gL < 

1 2 



TV 



\n{Ec/A) 



< 9L <1 , gR< fjL 



A 

'e^ 



9L 



Ec 



gL tan 



Ec\ 
A ) 



exp 



TT 




exp 
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